Fully permanent magnet atom chip for Bose-Einstein condensation 



OO 

o 
o 



(N 



5h 



i 

C 

o 
o 



> 

\o 
in 

(N 

in 
o 

o 



X 



T. Fernholz, 1 ^ R. Gerritsma/'Q S. Whitlock, 1 - 2 I. Barb, 1 and R. J. C. Spreeuw 1 

1 Van der Waals-Zeeman Institute, University of Amsterdam, 
Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 
2 ARC Centre of Excellence for Quantum- Atom Optics and 
Centre for Atom Optics and Ultrafast Spectroscopy, 
Swinburne University of Technology, Hawthorn, Victoria 3122, Australia 
(Dated: March 24, 2008) 

We describe a proof-of-principle experiment on a fully permanent magnet atom chip. We study 
ultracold atoms and produce a Bose-Einstein condensate (BEC). The magnetic trap is loaded ef- 
ficiently by adiabatic transport of a magnetic trap via the application of uniform external fields. 
Radio frequency spectroscopy is used for in-trap analysis and to determine the temperature of the 
atomic cloud. The formation of a Bose-Einstein condensate is observed in time-of-flight images and 
as a narrow peak appearing in the radio frequency spectrum. 

PACS numbers: 33.20.Bx, 37.10.Gh, 81.16.Ta 



I. INTRODUCTION 



Miniature patterns of magnetic field sources integrated 
on planar substrates, 'atom chips', are used to manip- 
ulate magnetically trapped clouds of ultracold neutral 
atoms PHI, IS IU- Atom chips based on current-carrying 
wires are relatively simple and versatile tools, readily 
used to produce Bose-Einstein condensates (BEC) 0, [g] . 
Subsequently a wide range of experiments have been per- 
formed; for example, BECs can be precisely positioned 
with atomic conveyer-belt potentials 0, 0] and phase- 
coherent splitting and interference of condensates has 
been achieved using double- well potentials [1, [§] . 

Atom chips can also incorporate permanent magnetic 
materials [l(| HH, [H, [H, [H[ , with significant additional 
advantages over current-carrying wires alone. Specifi- 
cally, permanent magnet atom chips do not suffer from 
current noise which is partly responsible for trap loss and 
heating of atoms near conducting surfaces [10j. They are 
also well suited for tightly confining trapping potentials 
without the technical problems associated with high cur- 
rent densities. Furthermore, patterned permanent mag- 
netic materials provide greater design freedom, allowing 
complex magnetic potentials to be realized such as ring- 
shaped wave guid es or large two-dimensional arrays of mi- 
crotraps EMUS- 

One obvious drawback of permanent magnets however 
is the relatively limited degree of dynamic control avail- 
able over the associated trapping potentials. The per- 
manent magnetic field is a particular source of difficulty 
when loading the atom chip. As such, most permanent 
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magnetic structures rely on current carrying wires to load 
and to complete the permanent magnetic trapping poten- 
tial. Analysis is another complication, because standard 
techniques such as free ballistic expansion are impossi- 
ble in the permanent magnetic field. These basic chal- 
lenges need to be addressed before more sophisticated 
atom chip designs incorporating permanent magnets can 
be realized. 

In this paper we describe the loading of ultracold atoms 
and evaporative cooling used to produce a BEC in a fully 
permanent magnetic microtrap, without additional ex- 
ternal control fields. Although the permanent magnetic 
chip is in principle " self-biasing" we use a small field de- 
rived from an externally installed permanent magnet to 
tune the Ioffe field of the trap. To load the permanent 
magnet trap we adiabatically transfer the atoms to the 
chip surface from a quadrupole magnetic trap by apply- 
ing a combination of bias magnetic fields. During load- 
ing we apply a radio frequency field to evaporatively cool 
the cloud. Analysis of the cloud temperature and spec- 
tral distribution is then performed using radio frequency 
spectroscopy QJ, [H Hfl SJ [H, H • This technique is 
also used to observe the BEC transition in-situ without 
expansion from the trap. It is also possible to eject the 
cloud from the self-biased trap by rapidly switching on 
uniform magnetic fields to observe the BEC transition 
using conventional absorption imaging. 



II. SETUP 

Our atom chip consists of a F-shaped permanent mag- 
net structure (Fig. [T]) , designed to produce a self-biased 
Ioffe-Pritchard microtrap 180 fim from the chip sur- 
face [H H3]. It is cut from 40 //m thick FcPt foil by 
computer numerical controlled (CNC) spark erosion and 
is glued to an aluminium-coated glass substrate. The 
magnetization is oriented in-plane along the y-direction 
and was measured using superconducting interference de- 
vices (SQUID) magnetometry as 430 kA/m. Assuming 
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FIG. 1: Permanent magnetic structure for producing a Ioffe 
trap 180 /J,m above the surface. The structure has a thickness 
of 40 /im and an in- plane magnetization of 430 kA/m. 



a uniform magnetization, we calculate an anisotropic, 
harmonic trapping potential with trapping frequencies of 
27rxll kHz and 27rx30 Hz in the radial and axial direc- 
tions respectively. The calculated magnetic field at the 
trap minimum (Ioffe field) is 1 G, in the x-direction. 

It was found experimentally that uncompensated stray 
fields, such as the earth magnetic field in combination 
with field corrugation caused by the coarseness of the 
magnetic structure, added field components in the x- 
direction on the order of 2 Gauss. The design aimed at a 
rather low Ioffe-field to achieve tight radial confinement, 
making it necessary to tune the Ioffe-field to a higher 
value to avoid field zeros in the trap. This was done with 
an externally installed permanent magnet, thus preserv- 
ing the permanent magnetic nature of the trap. The per- 
manent magnet was placed on the optical table (« 25 cm 
below the chip) , adding a homogeneous field of w 4 G in 
the x-direction. The trap frequencies in the lowest min- 
imum were then measured to be 27rx5.2 kHz radial and 
27rx 128 Hz axial with a field minimum around 2.7 G. Due 
to the field corrugation, the axial confinement is signifi- 
cantly larger than expected. We attribute the field corru- 
gation mainly to the rather unsophisticated production 
of the chip and feel that a more advanced atom chip (for 
instance, produced with patterning techniques based on 
optical lithography) could be fully self-biased. We note 
that in fact it was possible to load a secondary axial mini- 
mum of the trap and perform evaporative cooling without 
the permanent magnet installed, but less efficiently. 

The atom chip is mounted upside down in a quartz-cell 
vacuum chamber (4x4x7 cm 3 ), which is surrounded 
by three orthogonal coil pairs that can produce magnetic 
fields > 100 G with any desired polarity. The coils are 
used to load the atoms into the permanent magnet trap 
as described below. The loading procedure ends by ramp- 
ing the currents to zero and short-circuiting the coils with 
relays. The latter significantly reduces the residual mag- 
netic field noise. Additional sources of magnetic field 
noise, such as the current supplies themselves, were iden- 
tified and placed > 5 m from the experiment. Residual 
background field noise at 50 Hz was further reduced by 
feeding a phase and amplitude adjusted 50 Hz signal de- 
rived from the mains to a large coil pair placed around the 
experiment. Our efforts reduced the root-mean-square 
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FIG. 2: Homogeneous magnetic field components applied dur- 
ing the loading process. Position I corresponds to the initial 
coil based quadrupole trap. The magnetic field zero is moved 
close to the chip and the anti-Helmholtz coils are subsequently 
turned off near position II, turning the trap into a Ioffe trap. 
This trap is then moved towards the final location, replacing 
the initial self-biased trap. All fields are turned off at position 
VI. 

magnetic field noise in the x direction from more than 
20 mG to below 0.5 mG. Radio frequency (rf) fields used 
for forced evaporative cooling and for rf spectroscopy of 
the atom cloud are produced by a computer controlled 
10-bit direct digital synthesis rf generator [25[ (AD9959) 
with four independent outputs. We use two outputs con- 
nected to separate antennas for evaporative cooling and 
rf spectroscopy. 



III. LOADING PROCEDURE 

In this section we describe the loading process of the 
permanent magnetic trap using external uniform bias 
fields. The possibilities and restrictions for moving mag- 
netic traps by this method have been described in detail 
in [26|. Efficiently loading atoms into the trap proves 
to be nontrivial since at the end of the loading proce- 
dure all external fields should be turned off. The main 
difficulty arises from the already present field minimum, 
located close to the chip, containing no atoms at the start 
of the loading process. This minimum can act as a leak 
when the trap containing the atoms is brought too close 
to it. We have optimized a particular loading trajectory 
that avoids this secondary minimum as much as possi- 
ble. Another important feature is that the magnetic field 
minimum is kept above zero during the last part of the 
loading process, making it possible to evaporatively cool 
during the transfer without significant loss due to Majo- 
rana spin-flip transitions. 

We collect approximately 7 x 10 7 87 Rb atoms in a mir- 
ror magneto-optical trap (MOT). The MOT is made with 
the help of external field coils, at a distance of 8 mm 
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FIG. 3: (Color online) Simulated contours of magnetic field 
strength [corresponding to Fig. 011- VI)] together with their 
x- and ^-projections (red=8 G, green=16 G, blue=24 G) dur- 
ing the loading process. The filled trap on the left-hand side is 
compressed and moved towards the chip while the empty trap 
is decompressed, moves away from the chip, and disappears 
when the field is turned to zero (VI). 



from the surface where the field of the chip is negli- 
gible. The atomic cloud is then compressed (CMOT) 
and optically pumped to the magnetically trappable 
\F = 2, nip = 2) state. After switching all laser beams 
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FIG. 4: (Color online) Spectrally colored absorption images 
taken during loading without evaporative cooling. The field 
of view is 3.6 x 1.7 mm 2 . Every image is composed of several 
absorption images taken at 22 different probe detunings, rang- 
ing from to 60 MHz. The absorption images are overlapped 
after mapping them onto a color code where red corresponds 
to G Zeeman shift and blue to ~ 30 G Zeeman shift, (a)-(d) 
correspond to Fig. [3] (II, IV- VI). (d) corresponds to the loaded 
final trap with all externally controlled fields turned off. 

off, the cloud is trapped 8 mm below the chip surface us- 
ing the quadrupole magnetic field produced by one pair 
of coils operated in anti-Hclmholtz configuration. Ad- 
ditional uniform fields can bring this trap closer to the 
chip. In close proximity to the chip surface the field gra- 
dient produced by the F-shaped magnet is large enough 
to hold the atoms against gravity and the anti-Helmholtz 
coils can be turned off. At this point two magnetic field 
minima exist. A weak trap containing the atoms is pro- 
duced by the magnet in combination with a uniform bias 
field and is located far from the chip surface. A tighter 
trap close to the chip surface also originates from the 
magnet and is at this point still vacant. 

An efficient method to transfer atoms to the permanent 
magnetic trap is to move the atoms to the location of the 
final trap while avoiding the secondary trap. We use a 
time ordered sequence of uniform magnetic fields in the 
y- and z-directions as shown in Fig. O The merging 
of the two minima is prevented through an additional 
B z component of the applied field. The B z component 
is then inverted, as shown in Fig. [51 moving the trap 
containing the atoms towards the chip, while the empty 
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trap retreats. Finally, after ramping the fields off, the 
atom cloud is left in the permanent magnetic trap. 

Three-dimensional calculations of the magnetic field 
iso-surfaces during the loading are shown in Fig. [3)^11- 
VI). During this transfer the magnetic field gradient of 
the filled trap increases from ~ 60 G/cm to ~ 6800 G/cm 
which leads to significant heating of the atom cloud. Al- 
though the two field minima do not merge, they pass each 
other separated by a finite potential barrier of ~ 15 G 
which allows energetic atoms to spill over to the decom- 
pressing trap. This loss mechanism can take away up 
to 80% of the magnetically trapped atoms. However, 
it has been possible to load 2 x 10 6 atoms into the per- 
manent magnet trap. A series of absorption images at 
various detunings are taken during the loading in order 
to create a series of optical spectra of the trapped cloud 
and are shown in Fig. 21 The observed atom distribu- 
tions compare very well with the calculated field distri- 
butions in Fig. [3] The cloud temperature inferred from 
the optical spectrum in the self-biased trap is approxi- 
mately T ss 1 mK, yielding a peak phase space density of 
p = N(hu>/k B T) 3 ps 8 x 10~ 7 . Here, N is the atom num- 
ber, h — h/2ir is Planck's constant, Q is the geometrical 
mean trap frequency, and ks is the Boltzmann factor. 

To improve the starting conditions for the final evap- 
oration stage and to avoid loss due to the finite barrier 
height we apply forced radio frequency evaporative cool- 
ing already during the compression of the trap. Evap- 
orative cooling however requires a trap with non-zero 
minimum to prevent spin-flip loss. We choose a trajec- 
tory that always maintains a non-zero minimum during 
the final approach of the atomic cloud. It is non trivial 
that such a trajectory exists [26| . We have found that 
the trap turns into a Ioffe trap after the turn off of the 
anti-Helmholtz coils and remains a Ioffe trap during the 
rest of the loading process (steps II- VI in Fig. [3]). The 
preliminary evaporative cooling stage during the transfer 
consists of a 6 s linear rf ramp from 35 MHz to 8 MHz. 
At the end of the ramp all externally controlled fields 
are turned to zero and the cloud is trapped in the po- 
tential produced by the magnetic structures alone. After 
loading, the cloud contains approximately 10 5 atoms at a 
temperature of 60 pK yielding a peak phase space density 
of p i=s 2 x 10~ 4 . Although the number of atoms is rel- 
atively small, the high trap frequencies result in a mean 
collision rate of ~ 200 s _1 . A final rf evaporative cooling 
stage is performed with a linear rf sweep over 2.3 s from 
8 MHz down to ~1.9 MHz which lowers the trap depth 
from 600 to 10 pK. 



IV. RADIO FREQUENCY SPECTROSCOPY 

The most commonly used diagnostic technique for ul- 
tracold atoms is free ballistic expansion followed by ab- 
sorption imaging. In our case this is not possible due 
to the permanent character of the magnetization of our 
chip. Instead, we employ radio frequency spectroscopy 
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FIG. 5: Radio frequency spectra taken for various final evap- 
oration frequencies of 370, 170, and 110 kHz above the trap 
bottom fo, as indicated by the dashed vertical lines [(a)- (c), 
respectively]. The model used to describe the spectra is ob- 
tained by numerical integration of the Landau-Zener spin- 
flip probability assuming Bose-Einstein statistics (solid lines). 
The calculations consider the anisotropic trap geometry and 
give a Rabi frequency of Q ~ 2n x 0.8 kHz. The dashed lines 
indicate the prediction of the potential energy model using 
Eq. ([3]) with the same parameters. 



to characterize the ultracold atom cloud in the perma- 
nent magnetic trap. This allows us to directly probe the 
energy distribution of the confined cloud and to extract 
the cloud temperature without ballistic expansion. 

The cloud is probed by applying a radio fre- 
quency pulse of frequency / for a duration of 0.5 ms. 
Trapped atoms satisfying the resonant condition hf = 
gFPBB(x,y, z) may be coupled to untrapped magnetic 
states and are then immediately lost from the trap. Here, 
gp is the atomic Lande factor, ps is the Bohr magneton, 
and B is the local magnetic field modulus. We measure 
the number of lost atoms by applying a magnetic field of 
25 G in the y direction to push the atoms away from the 
chip surface and to weaken the trapping potential for de- 
tection. Resonant absorption imaging can then provide 
an accurate measure of the number of remaining atoms. 

In Fig. [5] radio frequency spectra are shown for various 
final evaporation frequencies. Above the BEC transition 
temperature T c , we can fit a theoretical model to the ra- 
dio frequency spectra to obtain the temperature of the 
cloud, whereas below T c radio frequency spectroscopy can 
be used to observe the BEC transition without free ex- 
pansion and to measure the chemical potential p. The 
following sections describe both techniques. 
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FIG. 6: Temperature inferred from the axial density distribu- 
tion (□) compared to temperature obtained from rf spectra 
(a) as a function of final evaporation radio frequency with 
respect to the trap bottom. The gray area gives the con- 
fidence interval for the temperature inferred from the axial 
density distribution, based on estimated systematic errors on 
the pixel size and lensing errors. The curve for T c is obtained 
from the atom number and the trap frequencies, and is cor- 
rected for the finite size of the cloud The dashed line at 
120 kHz corresponds to the observed onset of condensation. 



A. Thermometry 

If we assume that the atoms move through the res- 
onant surface at a constant speed perpendicular to the 
surface, we can describe the spin-flip transition probabil- 
ity by the Landau-Zener expression [U [28| . To obtain a 
model for the loss rate spectrum we integrate the spin- 
flip probability over the velocity distribution and over 
the resonant surface. We introduce the radio frequency 
above the trap bottom u> r { — 2ir(f — fa) and the thermal 
frequency w t h = ksT/h. In the case of a gas well above 
the condensation temperature T ^S> T c we find that the 
loss rate in an isotropic trap is well described by 



dN 
~dt 



1 



(1) 



where we introduced the resonant atomic spectral distri 
bution n WtS = dN/du> r f, using Boltzmann statistics: 
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As it can be seen from Eq. {!]), the loss rate becomes in- 
dependent of the trapping parameters and reflects the po- 
tential energy distribution [Eq. ^] for sufficiently small 
Rabi frequencies tt 2 <C uilo^/ ^/rr 



'. This arises from 
the fact that the out-coupling rate becomes independent 
of the momentum distribution under the above assump- 
tions. 

For temperatures approaching T c , Eq. |(5J) is no longer 
valid since it is based on Boltzmann statistics. For Bose- 
Einstein statistics we cannot evaluate the integrals to ob- 
tain an expression for dN/dt similar to Eq. (JXJ) . Instead 



we integrate the Landau-Zener probability numerically, 
including the trap anisotropy. To obtain the cloud tem- 
perature from our data, we fix /i using the measured atom 
number N by setting J n Utf duj T f = N. The results of this 
calculation are shown as solid lines in Fig. [5] 

We can alternatively assume that Eq. ((T|) also holds 
for Bose-Einstein statistics with the correct potential en- 
ergy distribution (dashed lines in Fig. [5]). For a non- 
interacting cloud with T > T c the resonant spectral dis- 
tribution is given by 



LU 3 



93/2 



(3) 



Here 173/2(2;) is the poly logarithm function with base 3/2, 
and fj, is the chemical potential. The inferred tempera- 
tures deviate by less than 15% from the numerically ob- 
tained results. 

As an alternative temperature measurement we also 
observe the axial density distribution of the cloud in the 
trap. To avoid problems in the absorption imaging due to 
the high optical density of the cloud we lift the magnetic 
field at the trap bottom by 50 G to weaken the radial 
confinement, and detune the probe laser by 18 MHz to 
the blue of the a + transition. By fitting the axial density 
distribution we obtain the temperature of the cloud. 

In Fig. [5] we compare the temperatures found from 
spectra and the inferred temperature from the axial den- 
sity distribution of the cloud for different final evapo- 
ration radio frequencies. Also shown is the confidence 
interval of the temperature based on the axial size. This 
interval is estimated for systematic errors such as lensing 
of the probe light by the dense cloud, pixel size calibra- 
tion errors, and finite imaging resolution. The plot shows 
excellent agreement between the two temperature mea- 
surements. 



B. Bose-Einstein condensation 

To observe the onset of Bose-Einstein condensation we 
increased the probing time to 10 ms and decreased the 
radio frequency amplitude by a factor of 3 to minimize 
power broadening. It took considerable effort in magnetic 
field noise reduction to obtain the resolution needed to 
observe the narrow condensation peak. In addition, it 
was necessary to maintain the same cycle for the applied 
magnetic fields throughout the experiment to minimize 
hysteresis effects in the magnetic material and surround- 
ing structures, which otherwise could lead to variations 
in the trap bottom of several 10 kHz. The observed spec- 
tra below T c are shown in Fig.[7J Lowering the final radio 
frequency results in a bimodal spectrum and eventually 
in a spectrum of a nearly pure condensate with a width 
of ~ 5 kHz. This agrees well with the chemical poten- 
tial derived from atom number and trap frequencies of 
n/rriFh = 4.3 kHz. Each spectrum is averaged 6 times 
and consists of about 70 data points, spread nonlinearly 
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FIG. 7: Spectra of (a) a cloud close to T c , (b) a partially con- 
densed cloud, and (c) a nearly pure condensate. The dashed 
lines indicate the position of the final radio frequency cooling 
"knife" with respect to the trap bottom (/o ps 1.85 MHz). 
The full width of the condensate spectrum is ~ 5 kHz, cor- 
responding to the chemical potential fi/rriFh. The bimodal 
fits, shown in light gray solid and dashed lines for the conden- 
sate and thermal part, respectively, are based on a heuristic 
model. 



over the spectral range to emphasize the bimodal nature 
and condensate peak. 



V. EXPANSION 

As an alternative way to identify the BEC transition, 
we ramp on, over 2 ms, a magnetic field of 50 G, which 
rapidly transfers the atom cloud 1 mm away from the 
chip surface. At this position all magnetic fields are sud- 
denly switched off and the acquired kinetic energy of the 
atoms is sufficient to overcome the remaining potential 
barrier. The cloud is accelerated further by gravity and 
the rapidly decaying field of the chip. After 3.5 ms an 
absorption image is taken to reveal a partially expanded 
atom cloud (Fig. [8]). As the final trap depth is reduced, 
the atom distribution becomes bimodal with a narrow 
central peak and a change in aspect ratio. Combining 
these absorption images with the radio frequency spec- 
tra we find the BEC transition at a trap depth of about 
11 /iK, with 7000 atoms at a temperature of T c ~ 1.6 fiK, 
which is close to the expected value of T c = 1.3 /iK. Con- 
tinued evaporation results in a nearly pure BEC consist- 
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FIG. 8: Absorption images taken after expansion of the 
cloud (10 times averaged). The atoms are first taken out 
of the self-biased trap by field ramps. Then the fields are 
suddenly turned off and the cloud falls and expands for 
3.5 ms. From left to right, the figures show a thermal cloud, a 
partly condensed cloud, and a nearly pure condensate of 1500 
atoms, corresponding to final rf frequencies of 1970, 1920, and 
1890 kHz. The lower graphs and insets show axial and radial 
cross sections, respectively, through the center of the cloud, 
with two-dimensional bimodal fits shown with dashed lines. 



ing of ~ 1500 atoms. 

Although this ejection method is useful for observing 
the phase transition to BEC it cannot be reliably used to 
obtain temperatures because of the non adiabatic nature 
of the expansion and the large magnetic field gradients 
during the expansion. It can however be used to obtain 
condensate fractions. 

We observe strong heating of the ultracold cloud 
(~ 6 /iK/s) during the final stage of evaporation. This 
is likely caused by three-body collisions as the esti- 
mated density at the trap center near T c is quite high, 
n(0) w 6 x 10 14 cm -3 . As a result, the lifetime of the 
condensate is ~ 15 ms without a rf shield and ~ 150 ms 
with a rf shield [3(J [M| . The quoted density corresponds 
to a three-body decay rate of ~ 3 s" 1 [3l|, [§3]. It is thus 
likely that three-body decay limits the final number of 
condensate atoms during the last evaporation stage. 

It is interesting to note that condensation is observed 
in the expanded cloud at higher evaporation frequencies 
than in the rf spectra. We suspect this is because of 
the long rf probing time necessary to obtain the desired 
resolution and signal-to-noise ratio, causing considerable 
heating during the radio frequency pulse. It is also possi- 
ble that the signal-to-noise ratio does not allow the obser- 
vation of small condensate fractions, because it requires 
the removal of a substantial fraction of atoms. 



VI. CONCLUSIONS 

In conclusion, we have presented solutions to the prob- 
lems of loading and analysis of an atomic cloud on a fully 
permanent magnetic chip. We have demonstrated an ef- 
ficient loading procedure relying on the adiabatic move- 
ment of Ioffe traps using uniform external fields, allow- 
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ing the application of forced evaporative cooling during 
loading. After loading, radio frequency spectroscopy is 
used to characterize the cloud. A model is derived for 
fitting rf spectra to obtain cloud temperatures. The fit- 
ted temperatures compare well to temperatures derived 
from fitting the in-trap axial density distribution. Noise 
reduction and radio frequency power reduction give the 
resolution needed to observe the onset of Bose-Einstein 
condensation in radio frequency spectra. This results in 
bimodal spectra and eventually in a spectrum of a nearly 
pure condensate. Alternatively, condensation is observed 
by ejecting the atomic cloud by switching external mag- 
netic fields. Together, these techniques have allowed us 
to analyze our atomic cloud in the permanent magnetic 
field. 
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